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Abstract 

We study the transport properties of an overdamped Brownian particle which is simultaneously 
in contact with two thermal baths. The first bath is modeled by an additive thermal noise at 
temperature Ta- The second bath is associated with a multiplicative thermal noise at temperature 
Tb- The analytical expressions for the particle velocity and diffusion constant are derived for this 
system, and the reliability or coherence of transport is analyzed by means of their ratio in terms of 
a dimensionless Peclet number. We find that the transport is not very coherent, though one can 
get significantly higher currents. 
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I. INTRODUCTION 



The second law of thermodynamics imphes that it is not possible to get useful work from 
a single heat bath at constant temperature. The perfect time-reversal symmetry (or detailed 
balance) which exists under equilibrium conditions causes the flux in either direction to be 
the same resulting in zero average flux. But advances in this area, especially in the context 
of biological systems, show that it is possible to induce a directed motion provided there 
exist some nonequilibrium fluctuations. These are usually provided by unbiased external 
input agents. Devices which exploit these fluctuations to generate a unidirectional flow have 
been termed as ratchets or Brownian motors. For a detailed discussion of these systems, 
see fl. 

The underlying thermal fluctuations are usually quantifled in terms of temperature. 
Though these fluctuations cannot be observed directly in macroscopic systems, the recent 
advances in nanotechnology and molecular biology allow one to devise molecular systems 
where they play a pivotal role. Thus, rather than trying to avoid these inherent fluctua- 
tions, they have been utilized for constructive purposes: Brownian motors serve as the best 
example of this. 

Many different classes of model systems exist in the literature, and detailed studies have 
been done with regard to transport and energetics in such systems. In our present work, 
we consider a model of Brownian motor where the system has a space-dependent diffusion 
coefficient D{x). The main emphasis in these systems is that the potential need not be 
ratchet-like. Examples of such systems arise in semiconductor and superlattice structures, 
molecular motor proteins moving along microtubules, etc. 

aa. 

A space-dependent diffusion coefficient D{x) could arise either due to a spatial variation 
of temperature or friction coefficient, or both. For a system having a space-dependent 
temperature [T(a;)] the system dissipates energy during its evolution differently at different 
places. This implies that the system is out of equilibrium. A similar effect arises in a 
medium with a space-dependent friction coefficient [ti{x)] in the presence of an external 



noise 



-8j. To be more precise, frictional inhomogeneity does not generate any current by 
itself - nonequilibrium fluctuations are crucial for transport in these systems. 

It has been shown by Biittiker [8,^ that, in the presence of D{x), the Boltzmann factor 
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exp[—V{x)/kBT] has to be generalized as exp[—ip{x)] where 

^(x) = -|'rfg^. (1) 

Here, v{x) = —fidV/dx is the drift velocity and /i is the mobility. For a homogeneous 
system, D{x) is given by Einstein's relation D = fiksT. The stability and dynamics of 
nonequilibrium system is greatly influenced by the above generalized potential. Here, one 
has to invoke the notion of global stability of states, as opposed to local stability which is 
valid for equilibrium states [7|. 

In the present work, we consider a Brownian particle in contact with two thermal baths, 
A and B. The underlying potential V{x) and the friction coefficient ri{x) are periodic in 
space, and separated by a phase difference other than and vr. We study the quality of 
transport in such systems. Here, the unidirectional current arises due to a combination of 
both r]{x) and the fluctuations present in the second thermal bath. 

Millonas [sl and Jayannavar jc] have analyzed the basic framework for particle transport 
in such systems. Subsequently, Chaudhuri et al. jsl have also studied this model with 
modifications in the bath parameters. However, there has been no systematic analysis of the 
dependence of the current on system parameters. Further, there is no clear understanding 
of the coherence of transport in these systems. In this paper, we undertake a detailed study 
of both these properties. 

The transport of a Brownian particle is always accompanied by a diffusive spread, and the 
quality of transport is affected by this diffusive spread. This property has been quantified 
via the ratio of current to the diffusion constant and is termed as the Peclet number, Pe. 
The Brownian particle takes a time t = L/v to traverse a distance L with a velocity v and 
the diffusive spread of the particle in the same time is given by ((Ag)^) = 2Dt. The criterion 
to have a reliable transport is that ((Ag)^) = 2Dt < L^. This implies that Pe= Lv/D > 2 
for coherent transport. The Peclet numbers for some of the models like flashing and rocking 



ratchets show low coherence of transport with Pe ~ 0.2 and Pe ~ 0.6 pi|, respectively. 
Experimental studies on molecular motors showed more reliable transport with Pe ranging 
from 2-6 JjJ . Our earlier work on inhomogeneous ratchets in the presence of a single heat 
bath but, subjected to an external parametric white noise fluctuation, showed a coherent 
transport with Pe of the order of 3 [l^. From this and other works [l3^, one concludes 
that system inhomogeneities may enhance the effectiveness of transport, though sensitively 
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dependent on physical parameters. 

The present paper is organised as follows. Section II gives the basic description of the 
systems, and also the analytical expressions for current and diffusion needed to evaluate the 
Peclet number. In Section III, we present detailed results for the dependence of current and 
quality of transport on system parameters. In Section IV, we conclude this paper with a 
brief summary. 



II. MODEL 



We begin with the equation of motion for a Brownian particle of mass m and position x{t) 
coupled simultaneously with (a) an additive thermal noise bath at temperature Ta, and (b) 
a multiplicative noise bath having a spatially- varying friction coefficient //(x) at temperature 
Tb- The particle moves in an underlyiiig periodic potential V{x). The equation of motion 
of the Brownian particle is given by p, l6| 

mx = -T{x)x - V'{x) + ^7i{x)iB{t) + Uit)- (2) 

The Gaussian white noises ^a(^) and ^_B(t) are independent, and obey the statistics: 

{Uit)) = 0, {Ut)U{t')) = 2VAkBTA5{t - 1'), (3) 
{^B{t)) = 0, {^BmB{t')) = 2TBkBTB6{t - t'). (4) 

The two noises together satisfy the fluctuation-dissipation theorem and r(x) = V a + ^ BVi.^) ■ 
Here (...) denotes the ensemble average. In the present work, we have chosen V{x) = 
Vo(l — cosx) and = r^oil — A cos(a; — 0)]. Here < A < 1, and we set A = 0.9 throughout 
for optimum values. The phase lag between V{x) and r]{x) brings in asymmetry in the 
dynamics of the system. This inturn leads to an unidirectional current even in the presence 
of spatially periodic (nonratchet-like) potential. 

When the friction term dominates inertia or on time-scales larger than the inverse fric- 
tion coefficient, one can consider the overdamped limit of the Langevin equation. This 
corresponds to the adiabatic elimination of the fast variable (velocity) from the equation of 
motion by putting p = x = 0. This approach is only applicable for a homogeneous system. 
For an inhomogeneous system, Sancho et al. [l^ have given the proper prescription for the 
elimination of fast variables. The resultant overdamped Langevin equation for Eq. ([2]) is 



given by [4| 



i [x) L^{xj i [xj i (a;) 

The prime on a function denotes differentiation with respect to x. Using the van Kampen 

Lemma 15j and Novikov's theorem 16|, we obtain the corresponding Fokker-Planck equation 



for the probabihty density P{x, t) of a particle to be at point x at time t as follows: 



dP{x,t) _ d fV'{ 



dt 



dx {T{x) r2(x) 



P + 



TaTa d 
T(x) dx 



P 
f(V) 



+ TbTb 



T(x) dx 



x] 



r x) 



-p 



(6) 



When the potential V{x) is positive and unbounded, the system evolves to a steady-state 
distribution Ps{x) characterised by zero current. Setting J{x,t) = in Eq. (JG]), we obtain 



Ps{x) = N exp[—ip{x)], 



where 



iIj^x) 



dq 



V'{q)T{q) 



+ 



B 



Ta 



TaTeV'^q) 



(7) 



(8) 



.TATA + TBTBviq) r(g) T^r^ + TBr£.?7(g) 
Here, iplx) is the dimensionless effective potential, and is the normalization constant 

For periodic functions V(x) and ri{x) with periodicity 2tt, a finite probability current is 

obtained. Following Risken 

current J as 



J 



ITj . one readily gets the expression for the total probability 
[1 -exp(-27r(5)] 



Jo^ dyexp [-ij{y)] P^^" dxexp [^{x)]/A{x) 



Here, 6 determines the direction of current in the system and is given by 



6 = ^p{x) — %Ij{x + 27r), 



and A{x) in Eq.( |9]) is given by 



A{x) 



TaTa + TBTBr]{x) 
T^x) 



(9) 



(10) 



(11) 



From Eqs. ([8]) - ffTOj) . we can see that the system is in equilibrium when Ta = Tb, leading 
to zero net current. But when Ta Tb, the system is rendered nonequilibrium and one can 
extract energy at the expense of increased entropy. It can also be seen that no current is 
possible when either or is absent. Again, if ri{x) is independent of x, there will be no 
net current flow in the system. When Ta — Tb changes sign, the current also changes sign 
but not the magnitude. 
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One can also obtain an analytical expression for the diffusion constant D by following 



19| as 



D 



r~'\dx/L)i4x) 



'go 



(12) 



where I+{x) and I-{x) are as follows: 



[X) 



IJx) 



A{x 
exp [— '?/'(x)] 



exp [ip{x)] 

x+L 



dy exp [-ij{y)] , 



x—L 



dy 



Aiy) 



exp [ipiy)] . 



(13) 
(14) 



Here, L represents the period of the potential and is taken as 27r. From the above equations, 
the velocity [v = 27rJ), diffusion constant {D) and the Peclet number (Pe) are studied as 
a function of different physical parameters. All the physical quantities are taken in dimen- 
sionless form. In particular, v and D are normalized by {Vq/tioL) and {Vq/tiq), respectively. 
Throughout our work we have set Vq, r/o and fc^ to be unity. Similarly, Ta,Tb and Ta,Tb 
are scaled with respect to Votiq and Vq respectively. We have used the Gauss-Kronord rules 



for numerical evaluations of the integrals involved 



2Q|. 



III. DETAILED RESULTS AND DISCUSSION 

In Fig. [H we examine the behaviour of the effective potential iplx) versus x, for two 
different values of (p. The parameter values are = 1, Tg = 0.1, = 0.1, A = 0.9 and 
= 0.9. For < < vr, the current is in the negative direction, while for tt < < 2tt, 
the current is in the positive direction. This can be directly inferred from the slope of the 
potentials. The tilt in the effective potential identically vanishes when = or tt. Hence it 
is expected that the unidirectional current does not arise for the case when the phase lag is 
or vr. In Fig. [H we also study the effect of changing Tb{= 0.1 or 0.9) on tplx) for a fixed 
= O.Svr value. We see that, as F^ increases, the tilt in effective potential increases. We 
also see from separate plots (not shown here) that the tilt is higher for smaller F^ values 
and also for higher temperature difference between the two baths. Hence, we set T4 = 1.0 
and Tb = 0.1 throughout this paper unless specifically mentioned otherwise in the captions. 

Next, let us study the behaviour of the probability current J, diffusion constant D and 
Peclet number Pe for various parameter values so as to extract an optimal set of parameters 
for coherent transport. 
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In Fig. |2]we plot J, D and Pe as a function of the phase difference for = 0.1, and 
Tb = 0.1, 0.9. The current and diffusion constant varies periodically with 0, as expected. 
We find that an increase in F^ causes a decrease in J and D, thereby causing an increase in 
the Peclet number. We can see from the plot of Pe vs. (f) that the transport is marginally 
coherent (~ 2) for a range of - values at the higher F^ - value. 

For a particular value of phase difference (0 = O.Svr), in Fig. [3] we plot J, D and Pe as a 
function of the temperature difference between the two baths, Ta — Tb- The other parameter 
values are Fa = 0.1 and F^ = 0.9. Here Tb is varied from to 2 for a fixed Ta = 1. As 
expected, the current is zero when the temperatures of the two baths are the same. We also 
observe that the direction of current reverses depending upon the temperatures of the two 
baths A and B. For Tb > T4, the current is in the negative direction and vice versa for 
Tb < Ta- Also, their magnitudes are different in both cases. 

From Fig. [3l we also see that while the current exhibits a peak at Ta — Tb — 0.5, 
D decreases with increasing Ta — Tb- On the other hand, Pe increases with increase in 
Ta — Tb, i.e., the transport goes towards the coherent regime. 

In Fig. m we plot J, D and Pe as a function of Tb, for F^ = 0.1 and F^ = 0.9. The 
other parameter values are specified in the figure caption. For F^ = 0.1 we see J ~ in 
the presence of a single heat bath (Tb — 0). As Tb is increased, the current starts to build 
up and shows a peaking behaviour. When the value of Tb — T4 there is no transport as 
expected. For a particular F^ value, this non-monotonous behaviour is seen only for lower 
values of F^- When the value of F^ is increased beyond a critical value, there is a non-zero 
current even when the temperature of bath B is zero. 

The plot of Pe vs. Tb in Fig. H] shows that, as F^ is decreased, the Pe value goes up and 
Pe — )■ 2 for small Tb when F^ = 0.1. (We have shown here only two representative F^ - 
values). We also see that the Pe value is more stable for lower F^ values, i.e., Pe ~ 2 up to 
some value of Tb, and then decreases with further increase in Tb- Thus, we conclude that 
Fa < 0.1 for higher transport coherence. 

In Fig. |5] we plot J, D and Pe vs. F^ for Fa = 0.1 and 0.9, = O.Svr, Ta = 1 and 
Tb = 0.1. In this case, for Fa = 0.1, the current shows a peak with F^, the diffusion 
constant decreases with increase in F^ and Pe — )■ 2 for F^ > 0.3. On increasing the Fa 
value, we see from separate plots (not shown here) that the values of J and D are reduced, 
but the value of Pe is almost the same, though Pe > 2 for much higher F^ values, say 
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= 5. Thus, Fig. [5] also leads to the conclusion that a lower Ta - value along with a higher 
- value aids transport coherence. 



IV. SUMMARY AND DISCUSSION 

Let us conclude this paper with a brief summary and discussion. We have studied the 
transport coherence of an overdamped Brownian particle in the presence of two thermal 
baths, namely A and B with different noise statistics. One of the baths, B, is characterized 
by the presence of a state-dependent friction coefficient. The system considered here is 
analogous to a simple model of a Maxwell-demon-type of engine, which extracts work out 
of the nonequilibrium state from thermal fluctuations by rectifying the internal fluctuations 
in the system 5, ^. The space-dependent friction coefficient in bath B is necessary for 
generating unidirectional current in the absence of a bias. The direction of current (but not 
its magnitude) changes sign when (T4 — Tb) changes sign. The current identically vanishes 
when the temperatures of both the baths are the same. Also, in the extreme high friction 
limit, the current vanishes as the particle cannot execute Brownian motion. 

We have systematically analyzed the behaviour of current, diffusion and transport coher- 
ence as a function of system parameters. The problem of coherence of transport has not 
been addressed in the context of these systems. Our work demonstrates that a combination 
of lower Ta - value and a higher Tb - value would lead to an optimum transport coherence. 
Our study of the Peclet number shows that, for most of the parameter space, Pe < 2. Thus, 
we conclude from the present study that the transport is not coherent though it is possible 
to get high current. 

Our study can be extended in several directions. For example, it would be interesting 
to study the efficiency of energy transduction in this Maxwell-demon-type heat engine by 
using the methods of stochastic energetics developed by Sekimoto 10]. The fact that the 
steady-state probability distribution Ps{x), Eq. ([7j) is a nonlocal function of V{x) and ri{x) 
itself can lead to much interesting physics and would form the basis of our future work. 
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FIG. 1: Plot of '(p{x) vs. x for the following parameter values, (a) Tb = 0.9 with (/) = O.Stt and I.Stt, 
(b) Tb = 0.1 and 0.9 with cf) = O.Stt. The other parameter values are = 1, Tb = 0.1 and 

= 0.1. 
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FIG. 4: Plot of J, D and Pe vs. Tb for = 0.1 and 0.9, F^ = 0.9, (f) = O.Stt, and Ta 
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FIG. 5: Plot of J, D and Pe vs. Tb for Ta = 0.1 and 0.9 with 
values are Ta = 1 and Tb = 0.1. 



O.Stt. The other parameter 
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